Finitely Correlated Generalized Spin Ladders 



A. K. Kolezhukfcl 

Institut fur Theoretische Physik, Universitdt Hannover, Appelstr. 2, 30167 Hannover, Germany 
Institute of Magnetism, National Academy of Sciences and Ministry of Education of Ukraine 
36(b) Vernadskii avenue, 252142 Kiev, Ukraine 

H.-J. Mikeskai 

Institut fur Theoretische Physik, Universitdt Hannover, Appelstr. 2, 30167 Hannover, Germany 
(March 13, 1998; printed out on February 1, 2008) 

We study two-leg S = | ladders with general isotropic exchange interactions between spins on 
neighboring rungs, whose ground state can be found exactly in a form of finitely correlated (matrix 
product) wave function. Two families of models admitting an exact solution are found: one yields 
translationally invariant ground states and the other describes spontaneously dimerized models with 
twofold degenerate ground state. Several known models with exact ground states (Majumdar-Ghosh 
and Shastry-Sutherland spin-i chains, Ameck-Kennedy-Lieb-Tasaki spin-1 chain, A-chain, Bose- 
Gayen ladder model) can be obtained as particular cases from the general solution of the first family, 
which includes also a set of models with only bilinear interactions. Those two families of models have 
nonzero intersection, which enables us to determine exactly the phase boundary of the second-order 
transition into the dimerized phase and to study the properties of this transition. The structure of 
elementary excitations in the dimerized phase is discussed on the basis of a variational ansatz. For a 
particular class of models, we present exact wave functions of the elementary excitations becoming 
gapless at second-order transition lines. We also propose a generalization of the Bose-Gayen model 
which has a rich phase diagram with all phase boundaries being exact. 
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I. INTRODUCTION 

In recent years, Heisenberg spin ladders have attracted 
considerable attention, particularly motivated by their 
peculiar properties of being intermediate systems be- 
tween dimensions one and two, as well as by the hope 
to get some insight into the physics of metal-oxide 
superconductors.tl It is now well established that "reg- 
ular" (i.e., with only longitudinal "leg" and transverse 
"rung" exchange couplings) S = \ isotropic spin ladders 
with even number of legs have a spin liquid ground state 
with short-range correlations and an energy gap, while 
odd-legged ladders have a quasi-long-range ordered gap- 
less ground state. 

On the other hand, "generalized" ladders includ- 
ing other couplings beyond the simplest case of rung 
and leg exchange are interesting toy models which 
can interpolate between a variety of systems, exhibit- 
ing remarkably rich behavior.cTEl Diagonal interactions 
cause frustration and can change the structure of the 
ground state;BB recently it has been also shown that 
biquadratic interactions are important since they tend 
to produce dimerization and may lead to phase transi- 
tion into a "non-HaldaueZj spin liquid state with absence 
of magnon excitations.Q'El In real magnetic systems bi- 
quadratic terms can arise due_to effective spin-spin inter- 
action mediated by phonons.cj 

In the present paper we study the class of generalized 
S = h ladder models which admit an exact solution for 
the ground state in terms of the so-called finitely corre- 
lated, or matrix product (MP) statesBtj The technique 



of MP states appears as a natural generalization of the 
well known valence bond solid (VBS) states ,Eil which al- 
lows one to construct rather complicated wave functions 
with given quantum numbersE3 and handle them easily. 
For standard Heisenberg spin models .-those MP states 
can serve as good trial wave functions ,B*E3 while for cer- 
tain "perturbations" of standard mpdfiki those wave func- 
tions become exact ground states-Effll] MP states have 
also proved to be a conve nient p.tool for variational study 
of elementary excitations Oii-§Q 

We use two MP ansiitze for the ground state wave func- 
tion having different translational symmetry properties: 
one is invariant under translation for one ladder rung, 
and for the other ansatz the elementary cell consists of 
two rungs. Respectively, we obtain two families of Hamil- 
tonians with exact ground states; one should emphasize 
that the Hamiltonian itself in both cases is translationally 
invariant, so that the second family describes models with 
spontaneous dimerization. (Those families were studied 
by us recentlyaQ in a somewhat more restrictive formu- 
lation.) The two families have non-empty intersection, 
which makes possible to determine exactly the boundary 
of a second-order phase transition from translationally 
invariant to dimerized phase. ■— . 

Thus we show that, in contrast to the common belief, Ej 
finitely correlated states can describe critical models. 
However, it turns out that in the dimerized phase the 
dimer order parameter Di = Si j • (Si^+i — Sii_i) and 
spin correlation functions behave in a rather peculiar 
way at the transition point r = r c (where r is cer- 
tain model parameter): (i) the spin correlation function 



1 



(5f5f + „) = Ase~ n /^ always decays exponentially, and 
the correlation length £ either does not exhibit any sin- 
gularities, or diverges as (r — r c )~ 2 , but in the latter 
case the prefactor As oc (r — r c ) vanishes exactly at the 
transition, thus preventing the appearance of the long- 
range spin order; (ii) the dimer order correlation function 
(DiDi+ n ) oc (r — r c ) 2 does not depend on the distance n 
and is just a constant smoothly vanishing at the transi- 
tion. 

For excitations in the spontaneously dimerized phase 
we propose a simple variational MP-type ansatz describ- 
ing the elementary excitation as a pair of solitons in dimer 
order; we show that the variational estimate for the gap 
goes to zero at the phase boundary. 

We show that for the certain class of transla- 
tionally invariant models having effectively a spin- 
1 Affleck-Kennedy-Lieb-Tasaki (AKLT)EI VBS ground 
state (formed by only triplet degrees of freedom at each 
rung) it is possible to write an MP-type wave function 
for the singlet excitation which is an exact eigenstate of 
the Hamiltonian and whose gap closes at two different 
second-order phase transition boundaries. For another 
class of models whose ground state is a simple product 
of singlet bonds along the ladder rungs, we give the con- 
dition of existence of the exact triplet excitation. 

We also consider in more detail a toy model being a 
matrix-product solvable "perturbation" of the "compos- 
ite spin" model first introduced by Bose and Gayencl and 
generalized later by Weihong et al.p in the same sense as 
the AKLT model is a perturbation of the standard 5 = 1 
Heisenberg spin chain. This model is remarkable because 
one can obtain exactly its full phase diagram, including 
several lines of phase transitions of the first and second 
order. 

The paper is organized as follows: in Sect. |ll] we explain 
the general procedure of findinga set of exact MP-type 
solutions, and sections [II and IV describe the two fam- 
ilies of solutions mentioned above; in Sect. |^ we also 
discuss the boundaries separating dimerized and non- 
dimerized phases and the structure of elementary exci- 
tations in the dimerized phase. Section [v] is devoted to 
the generalized Bose-Gayen model mentioned above and, 
finally, Sect. VI gives a brief summary. 



II. FINITELY CORRELATED EXACT GROUND 
STATES: THE CONSTRUCTION ROUTINE 

We start from a general form of the isotropic transla- 
tionally invariant spin ladder Hamiltonian with exchange 
interaction only between spins on plaquettes formed 
by neighboring rungs, H = X^^M+i: where the local 
Hamiltonian h is defined as 

hi,i+l = ^(JR + £ ii)Sl,i-S2 ) j + ~ e -R)Sl,'i+l'S2,i+l 



+ l // DD(Si ii -S 2 ,i+i)(S2,rSi i i+i) (1) 

+ ^ / Ri?X^l,i'^2,i)(Sl i i+l-S2 : i+l) — Eq , 

here the indices 1 and 2 distinguish lower and upper legs, 
and i labels rungs. The model is schematically repre- 
sented in Fig. [j] For periodic boundary conditions Jr is 
the coupling on the rungs, and Hamiltonians with differ- 
ent sr are physically indistinguishable; however, we will 
see later that it is necessary to introduce £r in order to 
include models with inequivalent legs. In addition, there 
are generally four different bilinear exchange couplings 
on the legs and diagonals and three biquadratic terms. 
Together with the constant term —Eq we have a total of 
ten parameters. For later convenience, wc will also use 
the following notation: 
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We look for the ground state wave function ^> in a 
form p£ a so-called finitely correlated^ or matrix product 
(MPjEi state, and we use the following two ansatze for 

K V = l[9i(u), ^r^n^n-lM'^M, (2) 



9i{u) 



u\s)i + \t )i -V2\t + i)i 
\/2|t-i)i u\s)i - \t )i 



(3) 



Here \s)i and \t^)i are the singlet and triplet states of the 
i-th rung, and «i, U2 are free parameters. Strictly speak- 
ing, <3>o in (|2|) is a matrix whose trace should be taken to 
get the wave function of a ladder with periodic boundary 
conditions. Three other linear independent wave func- 
tions tr(o-^^ ), t 1 = 0, ±1 describe states of a ladder 
with open ends having different behavior at the edges, 
which-jare similar to the edge states in the effective 5 = 1 
chain.oEj The ansatz ^) obeys rotational symmetry, i.e., 
tr($o) is a global singlet, and the states trfV^o) form 
a spin-1 triplet. The state ^ r Q lv is translationally in- 
variant under the translation for one ladder rung, and 
^>Q lm describes dimerized states with spontaneously bro- 
ken translational invariance (unless u\ = U2 when it ob- 
viously becomes identical to ^o^)- Translation of 'I'q 1 " 1 
for one rung leads to a different state with the same en- 
ergy. The construction (gl.was originally proposed as a 
variational wave function;EI later n the translationally in- 
variant form was used by usO to construct a class of 
exact ground states for the model (0) under somewhat 
restrictive assumptions £r = 0, Jl = J'l, Vll = Vdd, 
Vrr = 0, and we have reported several particular solu- 
tions of the type *|. m in Ref. 

One can observed that ansatze (||) describe several 
known examples of VBS-type states, e.g., *Q nv at u =J1 
yields the ground state of the effective AKLT chainEJ 
whose 5=1 spins are composed from pairs of 5 = ^ 
spins of the ladder rungs, 
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H = s » s «+i + ^(Sn • S, i+1 ) 2 , (4) 

n 

while for u = 1 and u = oo one obtains two degener- 
ate dimer ground states of the Majumdar- Ghosh (MG) 
modeljEj with singlets residing on the diagonals and on 
the rungs, respectively [the MG model results from ([j]) at 
Vrr,ll,dd = 0, J'p = 0, Jl = J'l, Jr. = Jd = 2Jl]- The 
wave function <J , g ilm at ui = — U2 = ±1 describes a state 
with checkerboard-type ordered singlet bonds along the 
ladder legs. Thus one may think of the construction (|^) 
as of interpolating between several types of VBS states 
mentioned above. 

The construction of exact ground states can be per- 
formed in the way outlined in Refs. [l0|,[l7| one has to 
require 'J'g to be a zero-energy ground state of the lo- 
cal Hamiltonian /i^i+i, which ensures that it is a ground 
state of the global Hamiltonian H. This yields the fol- 
lowing conditions: 

(i) for \Pg nv , the local Hamiltonian h\2 should anni- 
hilate all states being matrix elements of the product 
9i{u) ■ g 2 (u): 

hu{gi(u) ■ g2(u)} = 0, (5) 

while for 'J'q™ it is necessary that h\2 annihilates all 
states contained in the two products gi(u±) ■ 32(1*2) and 
9i ("2) • 32 (ui), 

h\2{g\{u\) ■ g2{u 2 )} = 0, 

hi2{gi(u2) ■ 32(1*1)} =0; (6a) 

(ii) all the other eigenstates of hyi have positive en- 
ergy. Then ^>q is the zero-energy ground state of H; if 
one drops the constant term — Eq in ([!]), the remaining 
Hamiltonian has the energy density Eq per rung. 

It is convenient to write the local Hamiltonian /ij^+i 
in terms of projectors on the states with fixed angular 
momentum \ipj k J) of the two-rung plaquette + 1). 
The complete set of the plaquette states contains one 
quintuplet (j — 2), three triplets (j = 1) and two singlets 
(j = 0), therefore a general form of h reads as: 

fj. k,l=l A 1 

k,l=l,2 

Here obviously Aj fe '^ = (Aj ^)* because of the hermitian 

property of h; ten independent A's are linearly related to 
the ten constants in the Hamiltonian of Eq. (|I]) . 

For each total momentum j, the complete set of the 
plaquette states can be divided into two subsets: rA local 

ground states \wf^) which are contained in the matrix 



products giQi+x, and n j local eigenstates ) which 

do not enter there; respectively the set of indices {k} for 
each j can be divided in two subsets which we will denote 
Gj for local ground states and £j for local eigenstates. 
The conditions (i) mean that the local Hamiltonian 

h should project only onto the states l^i^), and the 

multiplets \ipj 9 ^) have to be absent in Eq. (^|). This 
results in the following equations: 

Af'°=0, leGj,Vk (8) 

being essentially a system of linear equations in the 
Hamiltonian coupling constants J.., £.., V.., Eq [see Eq. 

(0)]. The only nonzero A's are then A^ for both k and I 
belonging to £j , so that hereafter we denote them simply 

as \f ,p ', p,p' = 1,.. .,n e 3 . 

The conditions (ii) require that all the eigenstates 
within the subspace determined by the set of states 
\ipj e ^} have positive energy, which yields the following 
inequalities: 

A<°> > , a = 1 . . . nj , (9) 

where Xj denotes the a-th eigenvalue of the matrix 

A, = [A^ p,) ], P y = l,...,nj 

If one or more of A^ is zero, this indicates a possibility 
for an additional degeneracy of the ground state (quan- 
tum phase transition). However, generally the conditions 
of stability imposed by the inequalities (|9|) are only suf- 
ficient but not necessary, and in order to conclude on 
the presence of a phase transition, one should be able to 
construct the other ground state explicitly. Later we will 
see that in many important cases this is indeed possible. 

III. MODELS WITH UNBROKEN 
TRANSLATIONAL SYMMETRY 

For 4 , nv , it is easy to check that generally the ma- 
trix product gi (u) ■ g2{u) contains only the following two 
multiplets which should become local ground states: 

(3 + u 4 y 1/2 (u 2 \ss) - V3\tt) Q0 ) = \Ao) . (10) 

(1 + u 2 )- 1 / 2 {^(M + K» - I**>im} = IV>S> ; 

(in fact, the second multiplet above enters only in the 
combination (3 + u 2 )^ 1 ^^ ), so that [V'iu) disappears 
from gi(u) ■ 32(1*) in the limit u — > 00, see below). The 
remaining, multiplets becoming local eigenstates can be 
chosen ascl 
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- (3 + u 4 )- 1 / 2 (^| SS > + u 2 \tt) 00 ) 



\4f) = J 7 n (Ks)-\st,)) 



V2 



(11) 



i<' 2) > = (i+- 2 r 1/2 (^(iv 
i4: ) > = i«> 2 p. 



i 



Here |ft)i M and |tt)oo denote respectively the triplet (j = 
1) and the singlet (j = 0) constructed from two triplet 
states on rungs i and i + 1. 

Equations (ph [in this case rig = rag = 1, nf = 1, raf = 
2, 7i| = 1, nf = 0], which impose the demand that the 
local Hamiltonian annihilates all the states (|l0|), can be 
rewritten in a form of relations between the Hamiltonian 
parameters as follows: 



AB 



1 



(l.i) 



(2,2h 



t(+) 



J 



(+) 

D 



5 X A{A + B) X 1 

A? — An — — 

12 12 4 

1^+5) 1 
12 2 12 4 V 1 



(A< M) +CA< 2 ' 2) ) 



l U.) _ ^(2,2) 



V LL : 



1 



Aq — Xi + CA^ 



: A 2 - 
3 

1 A(A + B) 

o A 2 H 5 A - A x 



(2,2) 



(12) 



CAi' 



(2,2) 



V - X A + ^i±^A ^A^ 1 ) A^ 2 > 



2Ai 



(1,2) 



wAi 



(1,2) 



2' 15 



<i n — 



wAi 



(2,2) 



the energy density per rung Eq being 

Eo = -J6 X2 -T6 X °- i^( A i 1,1)+A i 2 ' 2) ) 
Here we have used the following shorthand notation: 
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(3 + 1/2 ' 



B 



(3 + u 4 ) 1 / 2 ' 



C 



1 



1 



It is sufficient to consider only positive values of u, be- 
cause changing the sign u — > — u amounts just to inter- 
changing the ladder legs and thus does not bring in any 
new physics. One can see that the "unphysical" param- 
eter e R — — 2J^ '/it is just needed to allow the ladder 
legs to be inequivalent. The conditions (||) now take the 
form 

A o >0, A 2 >0, A^ 4) +A^ 2) >0, 
[X^KX^X^K (14) 

Equations (|l2| ) describe a hve-parametric family of 
spin ladder models (hve A's and u, but one of the A's 
is irrelevant since it just determines the energy scale) 
whose translationally invariant ground state can be writ- 
ten exactly in the form of a matrix product \&j, nv (g) . 



Though the limit u — > oo of Eqs. fll2] ) is formally well 
defined, the family ([l2]) does not contain all models hav- 
ing a product of singlet bonds on the rungs as their exact 
ground states, e.g., the model of Bose and Gayeno is miss- 
ing there. For u — oo obviously only the state \ss) enters 
the matrix product g\ ■ g2, and this case should be con- 
sidered separately. The resulting general solution for the 
u = oo case is 



J { l ] -J^ = {Vll-V dd )/A, 
the energy density per rung being 



(15) 



E Q = -3 J R /4 + 3{V LL + V DD + 3V RR )/16 . (16) 
The conditions of positivity (j^) result in the constraints 

IMj = Jr- 2J d +) + (V DD - V RR )/2 > 0, 
M2 = Jr + J [ d ] - (Vdd + 2V RR )/4 > 0, 



V? > 0, 

?(0 



1,2,3, 



(17) 



where are the eigenvalues of a symmetric 3x3 matrix 
with the matrix elements defined as 

ii _ Jr Vll 3V rr , 2 _ Jji Vdd 3Vrr 
Ail "2~2~4' Ml ~2~2 _ 4' 



Mi 



£r/2, Mi 



13 



-J 



(-) 
L ' 



/'l 



23 



J 



(") 



Mf = Jr ~ J L +) ~ (2 Vdd + 2V RR + V LL )/A . 



(18) 



One can check that the families ( |12| ) and (15) match 
each other if Ao = 0, so that it is natural to assume 
that vanishing Ao signals the first-order phase transition 
into the rung-dimer phase. Indeed, imagine we have a 
MP-solvable model with certain ground state and grad- 
ually decrease A to 0, then exactly at Ao = the rung- 
dimer state is degenerate with the "old" ground state, 
so that we have level crossing. The same reasoning can 
be applied to the transition into the completely polar- 
ized ferromagnetic phase which is driven by vanishing 
A 2 ; note that one cannot do the same for Ai because no 
uniform global wave function can be constructed from 
the local Hamiltonian eigenstates corresponding to those 
eigenvalues. However, later we will see that simultaneous 
vanishing of A^j 1 ' 1 ' and A^ can be associated with an- 
other phase transition, the second order phase transition 
into a spontaneously dimerized phase. In special cases, 
other eigenvalues may become releva nt for determining 
the phase transition points (see Sect. [II C and Sect. |V| 
below) . 

In the rest of this section we consider several particular 
solutions belonging to the general families derived above. 



A. Models with only bilinear exchange coupling 

If one demands that all biquadratic interactions in (|l|) 
are absent, the general family of models ( p"2| ) reduces to 
the following two "branches:" 
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1. Models with completely dimerized bonds, whose 
ground state is just a product of singlet dimers along 
the ladder diagonals, are obtained from (112) by setting 



i a (2,2) 

u = 1, X\ 
A 2 = 2 - y: 



1, Ai 



(i,i) 



(1,2) 



ty/2, A = 2y- 1, 



Jr = 2(1-2/), J L = l-y-x, J' L = l- y + x, 
J' D = 1, J D =0, £r = 2x, (19) 

the conditions Jl4| ) require that 1/2 < y < 2, a; 2 < j//2. 
It is easy to see that at x = this model describes an5 = 
i zigzag chain with alternated nearest-neighbor (NN) ex- 
change and frustrating next-nearest neighbor (NNN) in- 
teraction, with the alternation proportional to 2y — 1 and 
the ratio of NNN exchange constant to the smaller NN 
one equal to |. For ^ < y < 1 the ground state of 
this model was found by Shastry and SutherlandJHj and 
y — \ corresponds to the well-known Majumdar-Ghosh 
(MG) model;t3 at the MG point there is another, degen- 
erate ground state with singlet bonds on the rungs which 
is contained in the family (|l5|). At x = y — l-dl9|) de- 
scribes the so-called A-chain, or see-saw chainllHwhich, 
like the MG model, also has two degenerate ground states 
at y = ~. 

The interval 1 < y < 2 corresponds to partially ferro- 
magnetic interactions, and at y — 2 the model ([H]) ex- 
hibits a first-order phase transition into the completely 
polarized ferromagnetic state.S 

One can consider a more general model of a zigzag 
chain with two unequal NNN interactions a(l ± 7) and 
alternating NN couplings 1 and 0, < 1, as shown in Fig. 
||a. This model interpolates between the usual frustrated 
chain (7 = 0) and the A-chain (7 = 1). For fixed 7, the 
solution ([l9]) determines the line 

= 2a, max(-l, C-) < a < min(i, C+), (20) 

C± = {-l±(l + 8 7 2 ) 1 /2}/(47 2 ). 

in the (a0) plane where the ground state is a product of 
NN dimers on J = 1 bonds, see Fig. 0b. For | 7 | < 1 the 
minimal and maximal values of a in (EG) are always — 1 
and i, the point a = — 1, = — 2 lies on the boundary 
of the transition into the ferromagnetic (FM) phase, and 
a = 1/2, = 1 is the (generalized) Majumdar-Ghosh 
point. As we will see below, the situation for I7I > 1 is 
more complicated. 

It is worthwhile to make an attempt to estimate the 
boundaries of stability of the FM phase by using a stan- 
dard spin wave theory. Introducing two kinds of Holstein- 
Primakoff bosons for representing the spin operators on 
two different legs, and passing to the momentum space, 
one obtains in the lowest order the following spin wave 
Hamiltonian: 

H = ^{eifcflltflii: + £2ka\ k a 2 k + (®ka\ k a 2 k + h.c.)} , 

k 

1 + f 3 r, . u, f\ * P + e- lk 
£ik,2k = o a{l ± 7)(1 - cosfc), $ fc = . 



This Hamiltonian is trivially diagonalized, yielding the 
following expressions for the dispersion of two magnon 
branches: 



£A,B(k) = £ik, 



2 k 



(eife - £2k) + 2|$fc|2/b 



(21) 



(1 + 4) 

Zk = (£ik - £ 2fe )/2|$ fe | + (e lk - e 2fe ) 2 /4 1 ^ fc | 2 } 1/2 



It is easy to check that 
•i{2a + /3/(l +/?)}. Thus 



here a = sgn(e lfc - e 2k ) 
£ A (k) ~ Ck 2 at k -» 0, C -- 
£a represents the usual branch of ferromagnons, and the 
condition C = determines the point where this branch 
becomes unstable. The other branch eb represents "op- 
tical" ferromagnons, es(fc = 0) = — 1 — j3. For 7 = the 
branch eb is always gapped, but for large enough [7] it 
closes at k = n. Thus the region of stability of the FM 
phase is determined by the following two conditions: 



2a 



l + 2a 



H <7c(«,/3), 



(22) 



where the critical point j c is defined as a real positive 
root of the following fourth-order polynomial equation: 



64a 4 7* - [32a(2a + + 1) + 120 2 - 8/3 - 
- (0- 1) 2 (1 + 2a)(f3 + 2a) = 0. 



2„,2 



12K7 



(23) 



It is interesting to note that the curve in the phase 
space of the model determined by the first condition in 
( P2| ) at 7 = gives precisely the set of models with ex- 
act singlet ground states degenerate with the FM state 
which were studied in Ref. EO. Thus at least the first 



condition in (22), obtained from a simple spin wave cal- 
culation, remarkably coincides with the exact result. It 
should be mentioned that the exact singlet wave func- 
tions presented in Ref. |2^ exhibit double spiral order, and 
are very different from our MP wave functions, which in- 
dicates high degeneracy on the line of transition from the 
FM to singlet phase. 

The second condition in ( p2| ) becomes relevant only at 
I7I > 1, i.e., when the leg couplings have different signs; 
at | 7 | < 1 the boundary of FM phase is determined solely 
by the first condition, see Fig. 0b. 

One can also observe that the critical value of 7 de- 
termined by (p3) is different from that which can be ob- 
tained from our MP calculation; at a = —1, 8 = —2 it 
follows from ( pQ ) that |7| < 1, while from ( p3| ) one ob- 
tains I7I < 7 C ~ 1.7. As we discussed before, vanishing 
of any A's except Ao and X 2 cannot be associated with a 
phase transition, and generally the conditions (|l4[), be- 
ing sufficient but not necessary, tend to underestimate 
the stability region. 

2. Bilinear models with nontrivial ground states result 
from (112) when 



X 2 = u 2 (9-u 2 ), A = 2(3 + M 4 ), A^' 1} =5u 2 + 3, 
Af' 2) = (3 + M 2 )(l + « 2 ), A^ 2) ^x(l + u 2 )^ 2 , (24) 

which gives the Hamiltonian parameters as 
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J R = 6(u - 1), e R = 2x, 

J { L +) = u 2 {\ - u 2 ), = -ux, 

J D = -u(3 + u)(u- l) 2 , J' D = u(3-u)(m + 1) 2 , (25) 



2 (l+3u 2 ). From 



with the energy density per rung Eq - 
(H) it follows that x 2 < (3+u 2 )(5u 2 +3),"|m| < 3. At u 
3 the eigenvalue A2 vanishes, i.e., the system undergoes 
a phase transition into the completely polarized state, 
and at the same time one of the diagonal interactions 
vanishes, making the ladder equivalent to a frustrated 
chain;tl one can check that this solution belongs to the 
transition line from FM to the singlet phase determined 
by (^2|), and again the critical value of |7| obtained from 
the MP calculation, "f c = 1, is lower than one given by 
(§), 7c ~ 1.27. 

Another interesting limit is u — 1, when the rung cou- 
pling and one of the diagonal couplings are simultane- 
ously zero, and the ground state is a product of singlets 
sitting on the other diagonals. After a translation of the 
upper leg by one lattice site the model describes two spin- 
i chains, one antiferromagnetic and the other ferromag- 
netic, with exchange constants ± J, which are coupled via 
antiferromagnetic rung interaction J'; at J' J J > y2 the 
ground state becomes a trivial product of singlets along 
the rungs, a result which can be also obtained from (15). 



B. Multicritical models 

The Hamiltonian couplings in the general solution Jl2] ) 
become independent of the wave function parameter u if 



(2 2) ! 

Ao = X\ = X^'"' = 0, which defines a one-parametric 
family of models of the form 

J R = 2, J L = J' L = 1, J D = J' D = (1 + x)/{l - x), 

v 1/ 41-5z 41 + 5x 
Vrr = V ll = -- — , V DD = — , (26) 



(1,2) 



5 1 



5 1 



with a; = (3A^ /5A2) > 0. The ground state degeneracy 
of this model is very high: any wave function vE , o( u ) of 
the form (^]) with arbitrary parameter u is a ground state. 
One can easily calculate the overlap between two such 
wave functions having different values of u: 



(*o(«l)|*o(«2)) 



JV 



(1 + Ul u 2 ) 2 
(l + u 2 )(l + u 2 ) 



i.e., the two g.s. wave functions with different values of 
u are asymptotically orthogonal in the thermodynamic 
limit N — > 00 with the overlap vanishing exponentially 
with the increase of N. This means that the dimension 
of the basis of this subspace {^o(u)}, i.e., the number 
of mutually orthogonal ground states, is exponentially 
large in the thermodynamic limit. For x = this state- 
ment becomes trivial, because then the local Hamiltonian 
projects only onto the states with the total spin of two 



neighboring rungs equal to 2, and any sequence of sin- 
glets and triplets on the rungs which does not contain 
neighboring triplets is a ground state; in this case the 
SU(2) symmetry of the original Hamiltonian is sponta- 
neously broken because the ground states are not neces- 
sarily eigenstates of the total spin. Later we show that 
at x = I the model (|2fj) is gapless, by constructing the 
exact excitation wave function (see Sect. 0). 



C. AKLT-type models 



One can put u = in (|12j) and arrive at a large family 
of models having the AKLT ground state formed by ef- 
fective 5 = 1 spins made up of the ladder rungs. An 
example of such solution was presented in our earlier 
paper .0 The simplest model within this class can be ob- 
tained by requiring that the diagonal interactions van- 
ish together with the rung-rung biquadratic coupling, 

which is achieved by setting A 2 = 6/5, A^ 1,1 ^ = ^Xq — |, 

A 1 2 ' 2) = iA + § , A^ 1,2) = 0. The model then describes 
a ladder with ferromagnetic rungs and two additional bi- 
quadratic interactions: 

Jr. < -4/5, sr = 0, J L = J' L = 1, Jd = J' D = 0, 
V RR = 0, V ll = 12/5, V dd = -8/5. (28) 

It is worthwhile to mention that there exists a one- 
parametric family of modelaa connecting smoothly the 
Majumdar- Ghosh chain to the effective AKLT chain, 
thus proving that the Haldane chain and the dimerized 
chain are "in the same phase" in the sense that there 
is a continuous path in the phase space which links the 
two states and does not cross any singularities. However, 
as we will see below, this does not generally exclude the 
possibility to have a phase boundary between the Hal- 
dane (AKLT) and dimer states. A model exhibiting such 
a boundary will be considered later in Sect. |v|. 

The AKLT-type ladder models are especially remark- 
able because they admit constructing an exact eigenstate 
of the Hamiltonian being a singlet magnon which be- 
comes gapless under certain conditions. Let us consider 
the state 



W)s = tr{g^g^---gf l _ 1 g s n gf l+1 - 



■9n), 



(29) 



^ 1) ( 2 7) where gf = g^u = 0) is the ground state (AKLT) 



matrix and <?„ = g n (u = 00) = 1 • 



is the ma- 



trix describing a singlet bond on the n-th rung; a vi- 
sual VBS-representation of this state is shown in Fig. |j. 
This state is orthogonal to the ground state wave func- 
tion \&g nv = tr(<7i • • -g^). It can be straightforwardly 
checked that the action of the Hamiltonian on the state 
\n)s yields nothing but the states \n)s, \n±l)s'- 



H\n) s = {X^+X?< 2 >}\n) s 



(2,2)- 



+ -{\W-\^}{\n-l) s + \n + l) s }. (30) 
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[Here the ground state energy is already subtracted, 
i.e., the term — Eq is included in the local Hamiltonian 
h in (]]])]. Trivially constructing the exact eigenstate 
\k)s = J2n e1,kn \ n )s with total momentum fc, we obtain 
its energy as 



£s (k) = X^ + \™ + {X^> - Xl 1 ' 11 } cos(fe) . (31) 
The minimum of the dispersion is reached at k = 7T if 



v(2,2 



(2,2) 



A (i,i) < and at fc = o if \y-' L > > x\^>, and the ex 



,(2,2 



(1,1) 



,(2,2) 



citation becomes gapless when either Aj 1 ' 1 ^ or Xj*'* 1 van 



,(2,2) 



ishes. Later we will see that the transition at X\ 



(i,i) 



Ois 



into a spontaneously dimerized phase, and the transition 



at A 



(2,2) 



is specific for the AKLT-type models. 



Far from the phase boundaries one should expect that 
the singlet magnon branch will be high in energy, and the 
lowest excitation will be a^-usual triplet magnon being a 
soliton in the string order EM3 



D. Rung-dimer models with exact triplet excitations 

Inside the general family (|l^), (|l7j), describing models 
with rung-dimer ground state, it is worthwhile to con- 
sider a special subclass of models satisfying the additional 
condition 



J 



(-) 



(32) 



Those models are remarkable because it is a straightfor- 
ward exercise to check that the simplest excitation re- 
sulting from promoting one ladder rung from singlet to 
a triplet, 

\k) t =$> lfc » t (33) 

71 

\n) t = | s) 1 1 s) 2 ' ' ' \s) n -l\t)n\s)n+l ' ' ' \s)n, 

is an exact eigenfunction with the energy 
1 



st(k) = J R - -{Vll + V DD + W RR ) 



2(J« 



J^ +) )cosfc. 



(34) 



This exact excitation is of course not always the low- 
est one [e.g., for the model of A-chain, x = 1 — y in 
(|l9|), it is dispersionless and has a high energy equal to 
1], but in some cases the dispersion law (|I]) has gapless 
or almost gapless points and thus the triplet excitation 
defined above is relevant at least in some interval of mo- 
menta fc. One interesting realization is the ladder model 
with leg-leg biquadratic coupling, given by 



Jl = J' l = \v ll , Jr>U l >{), 



(35) 



with all the other couplings being zero. The dispersion 
( |34| ) in this case reads as e t (fc) = Jr — 2J^(1 — cosfc), 
and becomes gapless at fc = tt for — * \Jr. Another 
interesting model within this class will be considered later 
in Sect. [y|. 



IV. SPONTANEOUSLY DIMERIZED MODELS: 
PHASE BOUNDARIES AND ELEMENTARY 
EXCITATIONS 

For the dimerized matrix product ansatz ^/g lm the local 
ground states entering two matrix products g\(ui)- g 2 (u 2 ) 
and gi{u 2 ) ■ g 2 {u\) are 

\4o) = [3+K« 2 ) 2 ] _1/2 { W iU2|s S )-V3|tt)oo}, 

i^f 3 ) = [i + /72r 1/2 {(//2XK) + iv» - 



\Wil 2) ) = sgn(ui -u 2 ) 



V2 



\st L 



(36) 



where / = (ui + u 2 )/\/2, and the remaining multiplets 
can be chosen as 

\4o) = [3 + K U2 ) 2 r 1/2 {^l ss >+ u i"2|tt)oo}, 
- [2 + fr 1,2 {f\tt)i, + K) + M} 



(e)\ 
2ju / 



\tt) 



2/i- 



(37) 



Note that at u\ = u 2 the multiplet disappears 
from the set of local ground states and becomes an eigen- 
state, in accordance with Eqs. (pp|), (p"l|). Performing the 
same procedure as for undimerized case, we arrive at the 
following general family of solutions: 

1 ~~ 1 ~ 5 
J R = -ABX + -CAi + -A 3 , sr = 0, 
6 2 6 

Jl = J' l = -^A(A + B)X - IdXi + ^A 2 , 



jM = -±B(A + Bi\n--\ l ! -X,. 



1 

4' 



■5 

12' 



J 



(-) _ 2(ui + u 2 ) . 

(ui+u 2 ) 2 +4 u 



Vrr = ~ABX + CXx + ^A 2 , 
V LL = ^2(2 + B)X - CXi + h 2 , 
V DD = -B(A + B)X - Ai + i A 2 , 



(38) 



with the energy density per rung 

i^-lAo-lA.-lA,. (39) 

Here the quantities Xj > 0, j = 0, 1,2 denote the local 
eigenvalues corresponding to the multiplets and 
the factors A, B, C are defined as follows: 

uiu 2 + 3 ~ u\u 2 — 3 



.4 : 

C 



(m + tt 2 ) 2 -4 

(mi +m 2 ) 2 +4 



(40) 



Now we proceed to the study of general properties of the 
dimerized models derived above. 
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A. Phase boundaries 



B. Elementary excitations in dimerized phase 



Comparing the general solutions ( |l2| ) and ([38]), one can 
easily see that they match each other if 



Ml = U 2 



u, A 



.(1.2) 



(41) 



then X^' 2 ^, Ao, A2 in (|12|) s hould be identified respec- 



tively with Ai, Ao, A2 in (p8[). The original Hamiltonian 
(|l|) is translational invariant, thus the condition ([II]) de- 
fines a critical line in the space of the model parameters 
where the second-order quantum phase transition into 
spontaneously dimerized phase occurs. 

It should be mentioned that such a transition for lad- 
ders was first studied by Nersesyan and TsvelikQ within a 
field-theoretical approach in the approximation of weak 
inter-leg coupling. They found that leg-leg biquadratic 
exchange is a relevant perturbation, so that at sufficiently 
strong Vll the ladder can enter into a spontaneously 
dimerized non-Haldane spin liquid phase, with the ele- 
mentary excitations being soliton-antisoliton pairs. Later 
we have proposed explicitjexamples of solvable models ex- 
hibiting such properties ,0 and from the general solution 
( |38| ) it follows that other interactions may be responsible 
for this transition as well. (For example, it is possible 
to get a spontaneously dimerized ladder with Vll = 0). 
However, one can see that in MP-solvable ladder mod- 
els described by (^) the dimerization transition has cer- 
tain peculiarities: It is easy to calculate spin-spin and 
dimer-dimer correlation functions Cs{n) = (Sf A Sf i+n ) 
and C D {n) = (AA+n), A = S M • (S M+ i - S M _i) 
being the dimerization order parameter: 

C s (n) = {ul + ZY\Z 12 Z 21 ) n (& n . 2k + Z 2l 8 n<2 k+i), (42) 



C D (n) 



36(mi -m 2 ) 2 
(m 2 + 3) 2 (m 2 + 3) 2 ' 



Z a b = 



(u a -l)(u b + l) 



here a,b = 1,2. One can see that the dimer correla- 
tions exhibit long-range order vanishing for u\ — > m 2 
or M1.2 - ► 00, but remarkably there is no exponential 
tail present in the correlations. Spin-spin correlations 
always decay exponentially, and the correlation length 
£ = — 1/ \a.{Z\ 2 Zix) does not show any singular behavior 
at Mi — > M 2 , but diverges at u\$ — > 00; however, there 
is no long-range spin order at mi. 2 — * 00 since the ampli- 
tude of correlations in this limit vanishes. From the field- 
theoretical point of view, such peculiarities correspond to 
the "fine tuning" of the theory parameters which makes 
certain pre-exponential factors zero. 

It is easy to check that the general family of MP- 
solvable spontaneously dimerized models (|3^) includes 
one model of the multicritical typefl if u\ = —u 2 = u 
and Ao = 0, wave functions ^""(u) with any u have the 
same energy and are degenerate ground states. This so- 
lution can also be matched to the translational invariant 



family (jT|) at X\ 



(l.i) 



A = 0, m = 0. 



In order to discuss the structure of elementary exci-. 
tations in the dimerized phase, let us consider a modelEI 
defined by setting 

Mi = — M2 = m, Ao(u 4 + 10m 2 + 5) = 16/3, 

Ai = A (3m 4 + 14m 2 + 15)/8, A 2 = A (5m 4 + 18m 2 + 9)/8, 

which yields the Hamiltonian (|l]) with the following cou- 
plings: 



Vi 



DD 



-v. 



RR 



J, 



R 



8(m 2 -1)(m 2 + 3) 
3 m 4 + 10m 2 + 5 ' 



(43) 



4 5m 4 + 2m 2 + 9 1 

Fll = 3^ + 10^ + 5' Ji = J ^ = 1 ' J d = J d=0- 

The ground state energy per rung given by 



E Q = - 



7u 4 + 22m 2 
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4(m 4 + 10m 2 + 5) 



A "generic" example from this family is the model at 
u = ±1, when one gets purely biquadratic interchain in- 
teraction. It is remarkable first of all because of simplic- 
ity of its ground state, which is just a checkerboard-type 
product of singlet bonds along the ladder legs (see Fig. 
[|a); second, this model lies within the class, of Hamilto- 
nians considered by Nersesyan and Tsvelik.13 

Elementary excitations of the model (^3|) are pairs of 
solitons in the dimer order, and the scattering states of 
can be studied with the help o£ the following variational 
MP ansatz for a single solitonu 



C - "S^p i P( 2 ™+ 1 )|n\^ 



b>L = E e 



(44) 



N 



N< \t = Y[{92i-i(-u)g 2l (u)}g 2 '* +1 J| g 2l {-u)g 2i+1 (u), 

i—1 i— n+1 

gl = g(u) - Cs(-m), g\ tll = a"g(u) + Cfl(-«K ■ 

Here we for the moment assume that the ladder has 
2N + 1 rungs and periodic boundary conditions, so that 
the one-soliton wave function is well defined, u = 0,±1 
denotes the ^-projection of spin of the triplet excita- 
tion, and C = ±1 is the parity of a single soliton. The 
momenta arc defined in terms of the Brillouin zone of 
non-dimerized ladder, so that p € [0, w]. In case of the 
"generic" model with u = ±1 the states |n)f s can be 
visualized as singlet or triplet diagonal bonds separat- 
ing two "checkerboard-dimer" ground states, as shown 
in Fig. ||b. 

The lowest energy for the variational excitation (^) is 
always reached for the odd-parity state (£ = — 1), and the 
variational gap for the Haldane triplet is, always higher 
than the gap for a soliton-antisoliton pair.Q The energy of 
the elementary excitation (soliton pair) for the scattering 
states is given by 



8 



E(k,q) 



e s ,t 



[(k + q)/2] +e.,t[(k-q)/2] , 



where k and q are the total and relative momentum. 

The following expressions for the variational gaps can 
be obtained: 



Au 4 



{u 2 + 3) 2 



A 



(-) 



(u 2 + 3) 5 



(45) 



At u — > the odd-singlet soliton gap goes to zero, indicat- 
ing the second-order transition to the Haldane phase, and 
at u — > oo the odd-triplet soliton gap vanishes, signal- 
ing another second-order transition into the rung-dimcr 
phase. Thus, transition to the Haldane (AKLT) phase at 
u = is governed by a vanishing singlet-singlet gap, and 
the other transition to the rung-dimer phase at u — > oo 
is determined by the closing singlet-triplet gap; later in 
Sect. [v] we will see that the same property holds also for 
a different model, which suggests that it has a general 
character. 



V. GENERALIZED BOSE-GAYEN MODEL 

In this section we introduce a toy model which allows 
one to view simultaneously nearly all possibilities offered 
by the MP approach. Let us consider the model described 
by the Hamiltonian (|l|) with 

■h = J' L = 1, Jr = Vx, Jd = J'd = 2/2, Vrr = 0. (46) 

At j/2 — 1 and Vll = Vdd =t|0 this model was first 
introduced by Bose and Gayen,0 and its generalization 
to the case of arbitrary 2/2 W&s recently considered by 
Weihong, Kotov and Oitmaa.B 

The line 2/2 — 1 is quite peculiar because, as was shown 
by Xian,0 in this case the operator of the rung interac- 
tion commutes with the rest of the Hamiltonian, which 
enables one to classify the eigenstates by the total spin 
of each rung. At y\ > eq » 1.4 the exact ground state is 
just a product of singlet bonds along the rungs, and at 
yi < £o the ground state coincides with that of the effec- 
tive S = 1 Haldane chain whose S = 1 spins are formed 
by the pairs of S — h spins on the ladder rungs; Eq is 
exactly the ground state energy per spin of the S = 1 
Haldane chain, and at j/i = Eq a first-order phase tran- 
sition torn the Haldane phase to the rung-dimer phase 
occurs J3 When 1/2 7^ 1, this point of transition develops 
into a line as shown numerically by Weihong et al.El It 
is also worthwhile to mention that when the rung ex- 
change 7/1 is allowed to alternate, the systero-yxhibits a 
nice sequence of first-order phase transitions B3 

It is natural to try to "deform" the model by including 
biquadratic couplings Vll and Vdd in such a way that 
the Haldane phase becomes simply the AKLT phase, and 
then its ground state can be expressed through a matrix 
product wave function. It is a straightforward exercise to 
check that the general family (|lj) of MP-solvable models 



with translational invariant ground state reduces to the 
form Mq) under the following choice of parameters: 



Al 1 ' 1) = -^ 1 + |(4 y2 -l)>0, A^=0, 
2 5 



Ap ) = -> + i(4-2/ 2 )>0, u = 0, 
2 5 

A = -2/1 + ^(1 + 2/2) >0, A 2 = -(1 + 2/2) >0. 

Then the biquadratic couplings are expressed through 2/2 
as 



1 



(47) 



Vll = -(3 -2j/ 2 ) 



V, 



DD 



:(3lfe-2) 



(48) 



This model has the AKLT-type ground state with u = 
under the conditions imposed by the inequalities (|47 
with the ground state energy per rung 



E, 



AKLT 



1 

4 yi 



2/2J 



(49) 



From the general analysis presented in previous sec- 
tions we know that vanishing Ao corresponds to the phase 
transition into the rung-dimer phase, A2 = determines 
the transition into the completely polarized ferromag- 
netic phase, and A^ 1 ' 1 ' = gives the phase boundary 
of a transition into a spontaneously dimerized state, only 
the last transition being of the second order. The Hamil- 
tonian (Ej), ( f48| ) has a special symmetry (cf. Ref. ||): 
exchanging the two spins of every second rung is equiv- 
alent to exchanging the leg coupling 1 with the diagonal 
coupling 2/2 with simultaneous exchanging Vll and Vdd, 
so that there is a one-to-one correspondence between the 
eigenstates of the Hamiltonian at 2/2 > 1 and < 2/2 < 1, 
defined as 



m >->■ 2/1/2/2, 2/2 1-> 1/2/2, 2/2 > . 



(50) 



Exploiting this symmetry, one can obtain the "mirror" 
of the second-order transition line = in the other 

half-plane; this happens to be exactly the line defined 



by the equation 



(2,2) 



0. The line Aq = is invari- 



ant under the transformation (pOj). At the intersection 
(2 2) — 

of lines Ao = and A^ = we have the multicritical 
point (2/1 = 2, 2/2 = §) of the type . . . g(u) g{u) . . . corre- 
sponding to the model ( |2^ ) with x = | (see Sect. |III B| ) , 
and its "mirror" is obviously another multicritical point 
,2/2 = 1) of th e typ e . . . g(u) g{-u) g{u) g{-u) 



(fi . - , 

(see Ref. |7j and Sect. [IV A ), lying at the intersection of 



the lines A^ 1,1 ' = and An = 0. 



Since the model (|46|), fl48J ) belongs to the class of 
AKLT-type models considered in Sect. Ill, it is possible 
to write the exact wave function of the singlet magnon in 
the AKLT phase whose softening drives the second-order 
phase transitions at 2/2 = 3 + ig2/i and 2/2 = 4 — 
According to (|3l]), its dispersion law reads as 

E(k) = -2/1 + ^(1 + 2/2) + 2(1 - 2/2) cos(fe) , (51) 
5 
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so that the excitation becomes dispersionless at y% = 1. 



Near the lines A 



(1,1) 



0, A 



(2,2) 



those singlet ex- 



citation should be the lowest ones, and their condensa- 
tion determines the physics of second-order phase transi- 
tion. On the line of transition into the dimerized phase, 
2/2 = 5 + jqVi , the gap closes at k = it, and on the other 
line j/2 = 4— |j/i the gap closes at k = 0. 

The Hamiltonian (f46|), ( [48| ) satisfies also the special 
u = oo solution (p"5|), which means that it has a rung- 
dimer ground state within the region determined by the 
conditions (17), (jig); those conditions take the form 



^o = ^i 3) =yi-7(i + y 2 )>0, 



= -2/i + §(22/2-3) >0, 

Mi 2) = ^i + ^(2-3j/ 2 )>0, 
2 

M2 = V\ + -(I + 2/2) > 0, 
5 



1 



(52) 



and the ground state energy per rung for the rung-dimer 
phase is 



£o RD = -^1 + 1(1 + 2/2) 



(53) 



It should be remarked that the equations (|52|) indicate 
only the region of guaranteed stability of the rung-dimer 
phase, and generally they do not necessarily correspond 
to the real phase boundaries (recall that the conditions 
of stability imposed by the inequalities (|J) are sufficient 
but not necessary); for example, a naive application of 
the Eqs. j52|) to the original Bose-Gayen model would 
give j/ic = 2 as the boundary of the rung-dimer phase- 
while the correct region of stability, according to Xian,El 
is wider, y\ c 1.4. However, in the present case the lines 



/io = 0, juj = 0, Jlf = are true phase boundaries. The 
line /io = obviously determines the first-order phase 
transition because it precisely coincides with the Ao = 
line, so that we know exact ground states on both sides. 
Further, one can show that the lines ju^ = 0, Jl^ = 
determine second-order transitions. Indeed, the model 
(pi|), ( ^8[ ) satisfies the condition (|32|), and thus the wave 
function (33) is an exact excited eigenstate, with the dis- 
persion given by 

2 

s t (k) = yi - -(1 + y 2 ) + 2(1 - y 2 ) cos k . (54) 
5 

It is easy to see that the k = ir gap closes at the ju^ = 

(2) 

line, and the k — gap closes at the uX — line. 
Note that the symmetry transformation (|50|) leaves the 
dispersions ( pT| ) and (|54|) invariant, except for the change 
k 1 — ► 7r — k. 

Boundaries of the ferromagnetic phase can be obtained 
by requiring that A2 is the lowest eigenvalue of the local 
Hamiltonian (in this case only relative sign of A's mat- 
ters), which gives 



7(2) 



-2/1 + 4 < 2/2 < -1 ■ 

It is interesting to remark that the same conditions of sta- 
bility of the FM phase can be obtained from the standard 
lowest order spin-wave theory: one has two branches of 
ferromagnons with the dispersion laws 

S {+) (k) =-^(l+y 2 )(l-C0sfc), 

5 
6 

S (-\k) = --(1 + 2/2) ~ 2/1 + 2(1 - 2/2) cosfc , 
5 

and demanding the excitation energies to be positive one 
arrives at the same result. The transition at the line 
2/2 = —1 is of the first order (the entire branch e^ + '(k) 
of "usual" ferromagnons collapses), while for the line 
2/2 = 1 2/i + 4 the situation is different: the k = tt gap of 
the "optical" ferromagnons branch e^~'(k) closes, which 
suggests a second-order transition. 

Gathering up all that, one arrives at the complete 
phase diagram presented in Fig. ^|. It is interesting that 
the model (^), (f48|), despite its simplicity, has a rich 
phase diagram with five phases and three multicritical 
points. All transition lines are exact, and only inside the 
phases marked Dl and D2 the ground state is not known 
exactly. Existence of second-order transition boundaries 
to the AKLT and RD phases means that the symmetry of 
Dl, D2 is lower than that of the AKLT phase; this spon- 
taneously broken symmetry can be only a discrete sym- 
metry connected with the spatial parity of the system, 
which means that the ground state in Dl, D2 is twofold 
degenerate. The AKLT-D1 boundary matches our con- 
dition for the phase transition into the dimerized 
state with spontaneously broken translational symmetry, 
which indicates that Dl is dimerized. The properties of 
D2 are equivalent to those of Dl in the sense that the 
two phases are connected by the symmetry transforma- 
tion (bC 



VI. SUMMARY 

We have studied a class of generalized S — | lad- 
der models admitting exact solution for the ground state 
in terms of finitely correlated, or matrix product (MP) 
states. We use two different MP ansatze for the ground 
state wave function and respectively obtain two fami- 
lies of Hamiltonians with exact ground states, one being 
translationally invariant and the other one spontaneously 
dimerized. 

The two families have non-empty intersection, which 
enables us to determine the boundary of a second-order 
phase transition from translationally invariant to dimer- 
ized phase. We show that the behavior of the dimer order 
parameter and spin correlation functions at the transition 
is peculiar: The spin-spin correlation length £ can be ei- 
ther finite or diverge when approaching the transition 
point, but in the latter case there is no long-range spin 
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order since the preexponential factor vanishes exactly at 
the transition; the dimer order correlation function is just 
a distance-independent constant term vanishing at the 
transition. 

We show that for a certain class of translationally in- 
variant models (having effectively a spin-1 AKLT ground 
state formed by the triplet degrees of freedom at each 
rung) it is possible to write an exact wave function for the 
singlet branch of elementary excitations whose gap closes 
exactly at two second-order phase transition boundaries; 
one of those boundaries corresponds to the transition into 
dimerized phase, and the origin of the other one is differ- 
ent. 

For excitations in the spontaneously dimerized phase 
we propose a simple variational MP-type ansatz describ- 
ing the elementary excitation as a pair of solitons in dimer 
order; we show that variational estimate for the gap goes 
to zero at the phase boundary. 

As an illustration of the technique and the ideas of 
the present approach, we also propose a toy model be- 
ing an MP-solvable deformation of the "composite spin" 
model recently considered by Bose and GayenEI and later 
by Weihong et alB One can obtain exactly its complete 
phase diagram, including several lines of phase transi- 
tions of the first and second order. 
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FIG. 2. (a) The generalized zigzag chain with two un- 
equal next-nearest neighbor couplings; (b) its phase diagram. 
Thick dashed line j3 = F(a) = -2a/(l + 2a) and thin dashed 
lines j3 = / 7 (a) denote the boundaries of stability of the fer- 
romagnetic (FM) phase resulting from the first and the sec- 
ond conditions in (123), respectively; the actual boundary be- 
tween the FM phase and the singlet phase is determined by 
the two conditions /3 < F(a), (3 < f-,{a). The line j3 = 1, 
a < a c ~ 0.2411 is gapless, the rest of the singlet phase is 
gapped. For I7I < 1 the FM-singlet boundary is determined 
solely by the thick dashed line, and for I7I > 1 the FM phase 
starts to shrink with increasing I7I as indicated by the thin 
dashed lines. On the thick solid line (/3 = 2a) the ground 
state is a product of dimers along the J = 1 bonds. For 
I7I < 1 it is stable for — 1 < a < |, and (|, 1) is the gen- 
eralized Majumdar- Ghosh point; for I7I > 1, according to 
the criterion ( pp| ) obtained from the matrix product calcula- 
tion, the region of guaranteed stability of the dimer-product 
ground state on the line (3 = 2a shrinks with increasing [— ^ | as 
shown with the corresponding symbols (however, the MP cri- 
terion generally gives the conditions which are only sufficient 
but not necessary, so that the actual region of stability of the 
dimer-product ground state may be wider, see the discussion 
in the text). The point (— |,— |) at the thick dashed line 
(shown with an open circle) has a matrix-product state ^o nv 
of the type (Q) with u = 3, which is stable at least for I7I < 1. 
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n 

FIG. 3. Schematic representation of the exact excitation 
wave function ( p9| ) in the "AKLT-type" ladder model. Thick 
solid lines denote valence bonds connecting effective S = 1 
spins (indicated by ovals) which are formed by the triplet (t) 
degrees of freedom on the ladder rungs, and s denotes a singlet 
bond connecting two spins of the n-th rung. 
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FIG. 4. (a) Two degenerate ground states of the "generic" 
dimerized model (^) with u = 1; (b) the soliton states |n)^ 3 
used in Eq. (p4|), in the same special case u = 1. Thick solid 
lines indicate singlet bonds, and thick dashed lines can be 
either singlets or triplets. Arrows indicate the "direction" of 
the singlet bonds [i.e., |si^ 2 > = 2" 1/2 (| T1J.2) - | I1T2))]. 
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FIG. 5. Phase diagram of the generalized Bose-Gayen 
model @, (||). AKLT and FM denote the effective Af- 
fleck-Kennedy-Lieb-Tasaki (u = 0) and ferromagnetic phase, 
respectively, RD stays for the rung-dimer phase, and Dl, D2 
indicate phases with spontaneously broken symmetry where 
exact ground state is not known. The transitions AKLT-RD, 
AKLT-FM (shown with thick dashed lines) are of the first or- 
der, while all the other ones (shown with thick solid lines) are 
second-order transitions. Points Ml (|, |) and M2 (2, |) are 
multicritical with exponential degeneracy of the ground state, 
and M3 (—4, —1) is also a multicritical point but of different 
nature. The indicated eigenvalues of the local Hamiltonian 
vanish on the corresponding transition lines (on the FM-D1 
line A 2 = Aj M) ). 
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